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The problem of one-dimensional quantum wire along which a moving particle interacts 
with a linear array of N delta-function potentials is studied. Using a quantum waveguide 
approach, the transfer matrix is calculated to obtain the transmission probability of 
the particle. Results for arbitrary N and for specific regular arrays are presented. Some 
particular symmetries and invariances of the delta-function potential array for the N = 2 
00 . case are analyzed in detail. It is shown that perfect transmission can take place in a 

variety of situations. 
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1. Introduction 



The use of one-dimensional potential arrays is very frequent in several fields of 
physics. In solid state physics, the one-dimensional periodic arrays are used to model 
the lattice in a crystal, as a first approach. Several works use similar form ulatio n to 
demonstrate Anderson localization in different types of lattice disorder l l l ^ l ^K I 

A delta- function is useful to describe short-range potentials, like the interaction 
between the electrons and fixed ions in a lattice crystal. Thus, a periodic array 
of delta-function potentials is used in the Kronig- Penney model "I On the other 
hand, the delta-potential is also useful to describe impurities in solid state systems. 
Thus, the study of electron scattering by impurities in quant um wire s, using delta- 
potentials has been a subject of great interest in recent years EE1S19I j n so iir\ state 
quantum computation, finite ^-function potentials arrays are often us ed to desc ribe 
an instantaneous interaction between flying qubits and statics qubits I 10 | H | 12 | 

The problem of one-dimens ional delta-fu ncti on potential array has been studied 
in previus work j 1 3 | 1 4 | 1 5 | 1 6U71 

In reference they use a field-theoret ic approach 
to obtain the transfer matrix in terms of a propagator. In references I 14 | 15 | 16 LL7J 
the case of two delta-function potentials is studied using a quantum-mechanical 
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Fig. 1. (Color online) One dimensional quantum wire with N delta potentials of strength J n 
located at x = x n , n = 1, 2, .., N. 



approach. We use a similar aproach to develop a convenient way to deal with the 
problem, based in a transfer matrix methodology. 

We suppose a particle incident from the left on a linear array of N delta-function 
potentials. Each potential reflects and transmits a part of the particle wave function, 
which is taken to be a linear combination of an incoming and an outgoing plane 
wave in between two potentials. 

This paper is organized as follows. In section [2] the transfer matrix is obtained 
and some of its interesting properties are discussed. In section [3] results for specific 
regular arrays and arbitrary N are presented. In section|4]we analyze some particular 
symmetries and invariances of the delta- function potential array for N = 2, which 
has many interesting features. In that section we also analyze a condensed matter 
system: the transmission probability of one electron scattered by two impurities in 
a GaAs nanowirc. Finally, section [5] contains some concluding remarks. 



2. The transfer matrix 

We study the motion of a particle incident on N delta-function potentials located on 
a one-dimensional quantum wire along the £-axis. The Hamiltonian of the system 
is 



-2 N 

H= + Y\ J n S(x-x n ), (1) 
n=l 

where p and m are momentum and mass of the particle. The values x n and J n are 
the position and the strength of the n-th delta-function potential. The sign of J„ 
can be positive denoting a potential barrier or negative denoting a potential well. 

Figure [T] shows the array of the delta- function potentials. In this case an 
incoming particle from the left will be partially reflected and partially transmitted 
by each delta-function potential. As a result of these reflections and transmissions 
the particle wave function will be described by a combination of an incoming and an 
outgoing wave in between two potentials. The only exception is in the zone after the 
last (AT-th) potential. That is, in the region x > xn there will be only an outgoing 
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wave (e tkx ) since no particle is incident from the right. For each of the N + 1 zones 
in which the potentials divide the axis, the wave function is taken to be 

ip n (x) = A n e lkx + B n e~ lkx x n <x<x n+1 . (2) 

The wave amplitudes A n and B n are constant coefficients and k is the wave 
number for the particle energy e = h 2 k 2 /2m. This wave function is continuous at 
every point on the x-axis, however the delta potentials produce a discontinuity in 
the derivative of the wave function at the points x n , which is given by 

/ dip\ ( dip\ 2mJ n 

For the n-th potential J n S(x — x n ) at x = x n , the continuity of the wave function 
yields the relation 

A „lkx n , r> —ikx n a „ikx n , r> —ikx n j a\ 

A- n -ie + rs n _ie — A n e + ts n e . (4J 
Defining E n = e 2lkXn , Eq. [4] becomes 

A n E n + B n = A n -xE n + B n -i, (5) 
From the discontinuity in the derivative of the wave function at x — x n we obtain 



ik{A r 



„ikx n 
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') (6) 



or 



(1 + i\ n )A n E n - (1 - i\ n )B n = A„_i£„ - B n „ t , (7) 

in terms of the dimensionless parameter A„ = 2mJ n /h 2 k. Equations ([5]) and (J7J 
form a system of linear independent equations which can be expressed in matrix 
form as 



Bn I \ B 



>n-\ 



where 



2 V i\iE n 2 + i\ n 
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M„ is a transfer matrix which relates the coefficients of the wave function in the 
n-th zone with those in the (n— l)-th zone. As can be observed M n is a non- unitary 
matrix with determinant equal to 1 , it can be written in terms of the identity matrix 
I as 



M, = I ^L„ (10) 



where 



M^) = (_^5)- (ID 

Equation (|10[) divides M„ in a convenient way to separate the interaction. Further 
L„ only depends of the dimensionless parameter kx n . It is remarkable that this 
dependence is only in the off-diagonal elements of L„ , showing that the effect of the 
incoming wave over the outcoming one differs in only a phase. 

Now we can use Eq. ([8]) to relate the wave function coefficients of the outgoing 
wave function in the zone N, to the coefficients of the incident wave function to 
obtain 



(*:H (w >(*)' (i2) 

with Ai(N) = MjvMat_ 1 ...M2Mi. This formalism allow us to analyze the problem 
of a particle incident from the left on an array of N delta-function potentials. After 
the last potential, at x = x^, there is only an outgoing wave to the right, since we 
assume there is no particle (wave) incident from the right of the A^-th potential, 
that is Bn = 0. Since dct(M„) = 1 for all n and Bjy = 0, Eq. (TT2"j) gives 



An ~ jMikh, (13) 

The coefficient An is the transmission probability amplitude which means that the 
the probability of transmition is T = \An\ 2 /\Aq\ 2 . The incoming and outcoming flux 
have to be equal in each zone this means that \A n \ 2 — \B n \ 2 — \A n -i\ 2 — \B n -i\ 2 
for n = 1,2, N. 

We note here some interesting properties of the L„ which will be useful in the 
next sections. The anticommutator of the L-matrices is 



L„L m + L m L„ — (2 — E n E* n — E n E m )l. (14) 

Taking E n — E m , we note that \? n = which implies that M^ 1 = I + ^L n . Also, 
if E n = -E m then L„L m + L m L n = 41. Now, defining <f> mn = k(x m - x n ) Eq. (JHJ 
can be expressed as 
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L n L m + L m L n = 4 sin (0 mn )I. (15) 

Equation (|15[) is useful to simplify the multiple products of L's present in the transfer 
matrix. From Eq. (|T2j) . the general structure of the full transfer matrix Ai(N) in 
terms of A„L„ (n = 1, 2, .., N) is clear 



N N 

M(N) = I - - A ni L ni - - A ni L ni A„,L, 



2 ^/ ^ - '»ti — /ti ^ ^ ^ • -»ti — jii- v ri2 ^112 

ni—1 711,712 = 1 

ni>ri2 

AT 

2) / y A ni L ni A Jl2 L„ 2 ...A„ m L„ m + 

Hi ,n2 , . .,n m = l 
ri.i>n2>->n,T 1 

o ) A ril L ni A„ 2 L„ 2 ... A njv L njv . (16) 

If all A n are zero, that is no potentials are present, then obviously A4 = I. Equation 
(T17J)) simplifies enormously for special cases. For example, if all A„ = A for n = 
1,2, ..,N, A4(N) is polynomial of order N in the strength A. The information of 
the location of the potentials (x n ) appears in L„ through E n — e 2tkXn . Moreover, 
since products of L„ appears in Eq. I|16p it is through EnE^. Consequently, for 
specific regular arrays one can use Eq. (I15|) to simplify M. (N) when there are specific 
relations between the E n 's. 



3. Regular arrays 

Consider a regular array in which x n — x n -\ = d for n = 2,3, ...,N. The transfer 
matrix will be different for particles for different k. For a regular array the dimcn- 
sionless quantity kd plays the crucial role. We study a couple of specific examples 
for choice of kd. 

Case A. The simplest case is when kd = ir. This situation is known as resonant 
condition (RC) and is widely used to describe spin scattering EHHHHI. J n this case 
E n E^ +1 = 1, consequently all E n = E and L„ = L for n — 1, 2, 3, JV. Since 
L 2 = 0, Eq. 1161 reduces to simply 

N 

M(N)=l- l -Lj2*n. (17) 

n=l 

In this case A4 is linear and symmetric in the A„'s. The transmission probability 
(see Eq. (JT3J)) is 

1 4 
T = n = n ■ ( 18 ) 
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Moreover, if the potentials are both attractive (negative A) or repulsive (positive A) 
this can profoundly affect T, In fact T = 1 if the sum J3n=l ^« = ^' ^ ne individual 
values of the A„ do not matter as long as the sum is zero. 

Case B. Consider the case when kd = tt/2. In this case E n E* +1 = — 1. Conse- 
quently, E 2n+ i = Ei and E^ n = —E\ with n = 1, 2, .. . Correspondingly L2 n +i = Li 
and L 2n = Lj. From Eq. (JTSJ LiL] + L{Li = 41 and LiLi = L|lJ = 0, it can be 
shown that .M only contains linear terms in LiL} and LjLi for any given N. For 
N = 4 Ei = E 3 = E, E 2 = Ei = -E and L 4 = L 3 = L, L 2 = L 4 = L' . In this case 
A4(4) reduces to 

M(4) = (l-A3A 2 )I-|((A3+A 1 )L+(A4+A 2 )Lt)-i(A4A3+A4Ai+A 2 A 1 -A 3 A 2 )LtL 

+ ^(A 3 A 2 AiL + A 4 A 3 A 2 Lt) + iA 4 A 3 A 2 A 1 LtL. (19) 

Consequently 



i 1 

(M(4)) 2 2 = (1 - A 3 A 2 ) - -(A 3 + Ai + A 4 + A 2 ) - - (A 4 A 3 + A 4 A X + A 2 Ai - A 3 A 2 ) 

1 1 
+ 2( A 3A 2 Ai + A 4 A 3 A 2 ) + -(A 4 A 3 A 2 Ai). (20) 

We note that the case B can be reduced to case A by changing the wave number 
of the particle k — > 2k. Below we discuss the specific example of N = 2 in detail. 



4. Two delta-potential system 

We now focus on a N — 2 array. Although it is the simplest array, it presents 
interesting behavior like resonant tunneling. In this case 

K t Z ^2 T *Ai A2A1 

M(2) = I - — L 2 - — Li — L 2 L 4 , (21) 

and the transmission probability when the incoming flux has A = 1, (see Eq. (fl"3f 
)is 



\{M{2)) 22 \ 2 \(2 + iX 2 )(2 + i\i) + X 2 XiE 2 E* 1 \ 2 ' y ' 

Equation (|22|) shows that T depends only on the distance between potentials 
(x 2 — xi). Further, T is symmetric under interchange of Ai and A2 independently 
of the location of the potentials. The multiscattering present in the middle of the 
two potential can create a positive interference between the waves, and therefore, 
this can produce a resonant tunneling. Nevertheless, an extra phase is added to the 
wave function when it interacts with every potential increasing the complexity of the 
problem. To find the conditions when resonant tunneling can happen, we calculate 



October 15, 2010 0:17 WSPC/INSTRUCTION FILE text 



Transmission properties of the one- dimensional array of delta potentials 7 




Fig. 2. (Color online) T ex t as a function of Ai and A2 varying between —2 and 2. The white 
surface represents the case where n is an odd integer, and the gray surface represents the case 
where n is an even integer. Notice that for odd (even) n and Ai = A2 (Ai = — A2) we get perfect 
transmission (T = 1). 



the derivative of |(-M(2)) 2 2| 2 with respect to 021 = k{x2 — X\). The vanishing of the 
derivative gives the condition 



. 9 , 2(A! + A 2 ) 

tan202i = -; : — : — ■ (23) 

4 — A2A1 

Note that Eq. (|23|) is satisfied by 2(/) 2 i + mr, for integer n. This implies interest- 
ing behavior of T when the distance between the potentials is varied keeping the 
strengths Ai and A2 fixed. Keeping the first potential fixed at x = x\, let the po- 
sition of the second potential be at X2(n) — a; 2 + mr/2k, n = 0, ±1, ±2, .... In this 
case, using Eq. I23[ the extremal values of the transmission probability (T ext ) are 
given by 



4 + AiA 2 + (-l)«AiA 2 V( 2 + Ai) (2 + A 2 ) ' 

with A„ = \\ + 2. If n is an odd (even) integer, the value of T is maximum 
(minimum) when both Ai and A2 have the same sign. If Ai and A2 have opposite 
sign then T is maximum (minimum) for n even (odd) integer. The crucial point is 
the sign of the last term in the denominator of Eq. (|24|) . The behavior of T ext as 
a function of Ai and A2 is displayed graphically in Fig. [2] for even and odd integer. 
When the magnitude of the strength of both potentials is equal (Ai = A = ±A 2 ), 
Eq. ([24)) becomes to 
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Fig. 3. (Color online) Transmission probability T as a function of Ai when A2 = 4/Ai. The blue 
(solid) line shows the case when m is an odd integer and T is a maximum. The red (dotted) line 
shows the case when m is even and T is the minimum. 



lext- 8 + A 2 (4 + A 2)( 1± (_!)„)■ W 

For odd (even) n and Ai = A2 (Ai = — A2) we get perfect transmission (T = 1). 
Note that when we have an even n and Ai = — A2, then ^21 = rnn, and we obtain 
the case A described in section [3] 

The perfect transmition (T — 1) is present in the system where Ai = ±A2 due 
to the symmetry of the system. Notice that T — 1 with Ai = — A2 and 0i2 = nir 
independently of the values of the A's which implies that the extra phase added to 
the wave function in the scattering is equal to zero. In all other situations the extra 
phase depends on the values of the A's. 

For the particular case when A1A2 = 4, Eq. (f23f requires that cos202i = 0. In 
fact for 2021 = (2m + 1)tt/2, (m = 0, 1, 2...), E 2 E{ = exp[«(2m + 1)tt/2] and this is 
equal to (— l) m i. In this case, Eq. (f22|) reduces to 



12 + A2 + A2+4(-l)™(A 1 +A 2 )' v ; 

The condition A1A2 =4 requires both the potential strength parameters to have the 
same sign, even so one can obtain T — 1. In Fig. [3]T is plotted as a function of Ai 
when A2 = 4/Ai. One obtains T = 1 for Ai = A2 = 2 when m is an odd integer, and 
also for Ai = A2 = —2 when m is an even integer. This is really interesting because 
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Fig. 4. (Color online) Transmission probability T as a function of the electron energy in a two- 
potential array, when the effective mass of electron is 0.067 times the electron mass, the distance 
between potentials is X\ — X2 = 100 nm and Ai = A2 = 1 eV A(solid blue line), or Ai = — A2 = 1 
eV A(dotted red line). 

the potentials are both repulsive or attractive. The separation of the potentials is 
crucial in producing resonant tunneling. 

To set these ideas in a condensed matter system, we suppose that the two 5- 
function potentials represent two impurities fixed in a GaAs quantum wire, and we 
study the transmission of a ballistic electron, although the formulation in this paper 
can apply to any particle which interacts with dclta-function-like potentials. In this 
situation, the effective mass is 0.067 times the electron mass, and th e stren gths of 
interaction are J\ = J2 ~ 1 eV Aas reported for magnetic impurities EES. Figure 
[4] shows the transmission probability as a function of the electron energy if the 
distance between impurities is x± — X2 — 1000 A. 

The curves of T as a function of the electron energy in the cases when Ai = A2 = 
1 eV Aand Ai = — A2 = 1 eV Aare very similar, but the influence of the extra phase 
present in the first case (Ai = A2), produces a displacement in the values of the 
resonant energy. Notice that the peaks denoting the energies where there is perfect 
transmission (resonant energies) are more and more wide as the energy is increased, 
suggesting an already expected asymptotic behavior of T for perfect transmission 
for large values of energy. 

5. Conclusions 

In this paper we studied the problem of one-dimensional quantum wire along which 
a moving particle interacts with a linear array of N delta-function potentials using 
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a transfer matrix method. We showed that the transfer matrix (M„) for the n-th 
potential, can be expressed in terms of the unit matrix and a matrix L„ , which has 
interesting properties and useful anticommutation relations. The properties of the 
L„ matrix were used to calculate the transmission probability in the case of regular 
arrays. 

The case of the simplest array, namely N = 2, was considered in detail in sec- 
tion [4] This unexpectedly showed a variety of cases in which resonant conditions 
are possible. One obtained T = 1 even when Ai and A2 have the same sign. As 
shown above this is due to the interplay between the values of Ai, A2, the separa- 
tion (x2 — xi) between potentials and the wave number k of the scattered particle. 
We also show that T = 1 is only possible when Ai = ±A2 and for certain values 
of 021 = k(x2 — xi).The methodology developed here could be useful for the study 
of different potential arrays. We have obtained the transmitted wave function from 
all the information that defines the potential array (positions and strengths of the 
potentials). It could be also interesting to study the inverse situation: How much 
information about the potential array could be obtained from a scattering exper- 
iment? This could be useful in the implementation of a scattering-based quantum 
information system. 
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